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Abstract. By a refinement of the technique used by Johnson and Zumbrun to show stability under 
localized perturbations, we show that spectral stability implies nonlinear modulational stability of 
periodic traveling-wave solutions of reaction diffusion systems under small perturbations consisting 
of a nonlocalized modulation plus a localized perturbation. The main new ingredient is a detailed 
analysis of linear behavior under modulational data u' {x)ho{x), where u is the background profile 
and ho is the initial modulation. 



1. Introduction 

Stability and behavior of modulated periodic wave trains have received considerable recent atten- 
tion, both in the reaction diffusion and conservation law contexts; see, for example [SI, S2, DSSS, 
JZl, JZ2, JZN, BJNRZl, BJNRZ2, NRl, NR2] and references therein. The initial mathematical 
challenge of this problem is that the linearized equations, being periodic-coefficient, have purely 
essential spectrum when considered as problems on the whole line, so that there is no spectral gap 
between neutral and other modes, making difficult either the treatment of linearized behavior or 
the passage from linear to nonlinear estimates. 

This issue was overcome in the reaction diffusion context in the late 1990s by Schneider [SI, S2], 
resolving an at the time 30-year open problem. Using a method of "diffusive stability," Schneider 
combined diffusive- type linear estimates with renormalization techniques to show that, assuming 
"diffusive" spectral stability in the sense to be specified later of the background periodic wave, 
long-time behavior under localized perturbations is essentially described by a scalar heat equation 
in the phase variable, with the amplitude decaying more rapidly. 

More recently, Oh-Zumbrun and Johnson-Zumbrun [OZl, JZ2] in 2010 using rather different 
tools coming from viscous shock theory have resolved the corresponding problem in the conservation 
law case, showing that, again assuming diffusive spectral stability of the background periodic wave, 
behavior under localized perturbations is described, roughly, by a system of viscous conservation 
laws in the derivative of the phase and other modulation parameters, so that the phase decays 
at slower, errorfunction rate than in the reaction diffusion case. This description has since been 
sharpened in [JNRZ2] , showing that behavior consists of a modulation governed by the associated 
formal second-order Whitham averaged system plus a faster-decaying term. Further developments 
include the treatment of nonlinear stability of roll waves of the St. Venant equations [JZN, BJNRZl] 
and the resolution in [BJNRZ2] of the 35-year open problem of proving nonlinear stability of 
spectrally stable periodic Kuramoto-Sivashinsky waves. 

Applied to the reaction diffusion context [JZl], the approach of [JZ2] recovers and slightly extends 
the results of Schneider, yielding the same heat kernel rate of decay with respect to localized 
perturbations while removing the assumption of [SI, S2] that nearby periodic waves have constant 
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speed. However, one might hope to do better, based on the successful analysis in [JZ2], etc., in the 
conservation law case, of perturbations decaying at slower, errorfunction rate. 

As described clearly in [DSSS], formal (WKB) asymptotics suggest the same thought, yielding 
the description to lowest order of a scalar viscous conservation law 

(1.1) kt - uj{k)x = {d{k)kx)x 

for the wave number k = ipx, where frequency oj = tpt and ip is the phase. Here, Lo{k) = —c{k)k and 
wave speed c{k) arc determined through the nonlinear dispersion relation obtained from the mani- 
fold of periodic traveling-wave solutions u^{k{x — ct)), and d{k) through higher-order asymptotics. 
See [Se, 0Z2, NRl, NR2] for corresponding developments in the conservation law case. From (1.1) 
and the definition k = ipx, one might hope to treat initial modulations ■0|t=o = ho for which, not 
the phase ho, but its derivative dxfiQ is localized (L^), obtaining the same heat kernel rate of decay 
for k as found in the conservation law case [JZ2]. 

In this note, we show that this is indeed the case, proving by an adaptation of the methods of 
[JZl, JZ2] that diffusively spectrally stable periodic reaction diffusion waves are nonlinearly stable 
under initial perturbations consisting of a localized perturbation plus a nonlocalized modulation Hq, 
dxho G L^, in phase, with heat kernel rate of decay in the wave number and errorfunction decay in 
the phase. In a companion paper [JNRZl], we build on these basic estimates to establish, further, 
time-asymptotic behavior, validating description (1.1) in the long-time limit. Recall that (1.1) 
has been validated in [DSSS] in the related, asymptotically-large-bounded-time small-wavenumber 
limit and for data in uniformly-local Sobolev spaces"*^, in the somewhat different sense of showing 
that there exists a (5-family of sohitions of the full system S'^-close to a formal expansion^ in 5 
on intervals [0,T/5'^], where M and T are arbitrarily large and 5 is the characteristic size of the 
wavenumber of the modulation: that is, they build "prepared" data for bounded intervals rather 
than work with general ones globally in time. 

Consider a periodic traveling- wave solution of reaction diffusion equation ut = Uxx + fiu), or, 
equivalently, a standing-wave solution u{x, t) = u{x) of 

(1.2) ki^ut = kluxx + /(-") + k-^cux, 

where c is the speed of the original traveling wave, and wave number A;* is chosen so that 

(1.3) u{x + l)=u{x). 

For simplicity of notation, we will follow this convention throughout the paper; that is, all periodic 
functions arc assumed to be periodic of period one. 

We make the following standard genericity assumptions: 

(HI) / G C^(M), {K > 3). 

(H2) Up to translation, the set of 1-periodic solutions of (1.2) (with k replacing k^) in the vicinity 

oi u, k = ki,, forms a smooth 1-dimensional manifold {u{k, •)} = {u'^(-)}) c = c{k). 
Linearizing (1.2) about u yields the periodic coefficient equation 

(1.4) k^vt = k^Lv:={kldl + k^cdx + b)v, b{x) := df{u{x)). 
Introducting the family of Floquet operators 

(1.5) k^L^ := e-'^'^kM^^ = kl{dx + i^f + k^c{dx + ii) + 6, 



That is, following the notation of [DSSS], k = ipx € H^i with small norm := sup^^gg ||w||jys([2,_2,_|_i]) and 

s > M + 2 > 4. 

^But this expansion is only 5^-close to the approximate solution involving only the solution of the second order 
Whitham equation. 
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operating on periodic functions on [0, 1], determined by the defining relation 

(1.6) He'^'^f) = e'^'^iL^f) for / periodic, 

we define following [SI, S2] the standard diffusive spectral stability conditions: 

(Dl) A = is a simple eigenvalue of Lq. (Note that ^ = corresponds to co-periodic perturba- 
tions, hence Lq has always at least the translational, zero-eigenfunction u' .) 

(D2) (t(L^) C {A I 0?A < -0\i\^} for some constant 6' > 0. 

Then, we have the following Main Theorem, extending the results of [JZl] to nonlocalized per- 
turbations. 

Theorem 1.1. Let K>3. Assuming (H1)-(H2) and (D1)-(D2), let 

Eo ■= \\uo{- - ho{-)) - ■u(-)||Li(R)nff^(R) + ll^x^o||Li(M)nH^(M) 

be sufficiently small, for some choice of modulation Iiq. Then, there exists a global solution u{x,t) 
of (1.2) with initial data uq and a phase function ip{x, t) such that, for t > and 2 < p < oo, 

M- - i^i;t),t) - «(-)||l.(m) <m + ty'^^'-'M 

\\V.M-MLP(R)^Eo{l+t)--2('-yP\ 

and 

(1.8) ||?i(-,t) -u(-)||l-(r), ||V'(-,t)llL-{R) < ^0- 

In particular, u is nonlinearly (boundedly) stable in L°°(IR) with respect to initial perturbations 
vo = uq-u for which \\vo\\e := 'v^'i^d^ho£L^{R)c\HK(^) Eq is sufficiently small. 

Remark 1.2. It may seem more natural, and indeed is so, to introduce Hq and ip via 

v{x, 0) = uq{x) — u{x + /io(a:)), v{x, t) = u{x, t) — u{x + ip{x, t)). 

However, in doing so one introduces in the equation for v terms involving only il) and thus not 
decaying in time; see Lemma 5.1 below. For this reason we work instead with 

v{x, t) = u{x — i^{x, t),t) — u{x), 

that is, 

u{x, t) = u{Y{x, t)) + v{Y{x, t),t) 

where Y is such that 

Y{x, t) - ^(y (x, t),t) = X, Y{y - V(y, t),t)=y . 
Notice that we insure the existence of such a map Y by keeping, for any t, ||'i/'('5 OIIl°°(IR) bounded 
and ||V'a;(')i)||L°°(M) small. It should be stressed, however, that^ 

Y{x,t)=X + ^lj{x,t)+0{U{;t)\\L^\\M;mL^) 

SO that we are not so far from the natural (but inappropriate) approach. Furthermore, notice also 
that introducing the map Y above enables one to go back to the original unknown u{x,t) when 
desired. 

At a philosophical level, the main new observation here beyond what was shown in [JZl] is that 
linear perturbations that are "asymptotically modulational"^ behave as regards the linear solution 
operator S{t) as if they were asymptotically constant; in particular, dxS{t){u'ho) ~ S{t){u'dxho), 
where u is the background traveling wave (see Section 4). The analysis is based, at a technical level, 

^This follows from Y{X{y,t),t) - [X{y,t) + ^{X{y,t),t)] = il)(y,t)--4,{y-il){y,t),t). 

^In the sense that they converge at plus and minus spatial infinity to scalar multiples Hqu' of the derivative u' of 
the background traveling wave, corresponding at linear level to shifts u{x + /io(a;)). 
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on the decomposition of Floquet modes of u'Hq into the sum of Fourier modes of ho times the periodic 
function u'{x), and the separate analysis/recombination of each frequency range [— 7r+2j7r, 7r+2j7r]. 
We are unable to see such detail from the Green function description of [JZl, JZ2, JZN], nor from 
weighted energy estimates as in [SI, S2]. Indeed, this seems to be an instance where frequency 
domain techniques detect cancellation not readily apparent by spatial-domain techniques. 

In the process, we obtain an alternative proof of the basic estimates in [JZl], carried out entirely 
within the Bloch transform formulation without passing to a Green function description as in [JZl]. 
This elucidates somewhat the relation between the Bloch tranform based estimates of [SI, S2] and 
the Green function based estimates of [JZl, JZ2], in particular the role of the Green function 
(integral kernel) in obtaining the estimates of [JZl]; see Remark 3.2. 

Our analysis suggests that one might by an elaboration of the approach used here treat still more 
general perturbations converging asymptotically to any fixed periodic perturbations at ±00, not 
necessarily modulations, which would then time-exponentially relax to modulations following the 
dynamics of the corresponding periodic problem; see Remark 4.2.2. This would be an interesting 
direction for further investigation. 

We mention, finally, that the techniques introduced here are not limited to the reaction-diffusion 
case, but apply equally to the conservative case treated in [JZ2], yielding a comparable result 
of stability with respect to nonlocalized modulations [JNRZ3]. Indeed, this clarifies somewhat 
the relation between the reaction diffusion and conservation law case, revealing a continuum of 
models with common behavior lying between these two extremes. These issues will be reported on 
elsewhere. 

Remark 1.3. Without change, we may treat reaction-diffusion systems with a diagonal diffusion 
as considered in [DSSS]. More generally, all proofs go through whenever the linearized opera- 
tor generates a semigroup; in particular, our results apply to the (sectorial) Swift-Hohenberg 
equations treated for localized perturbations in [S2].^ It would be interesting to extend to the 
multi-dimensional case as in [Uj. 

Note: We have been informed by Bjorn Sandstede that similar results have been obtained by 
different means by him and collaborators [SSSU] using a nonlinear decomposition of phase and 
amplitude variables as in [DSSS]. 

2. Preliminaries 
Recall the Bloch solution formula for periodic-coefficient operators, 

(2.1) (Smix) := (e*^5)(^) = T e'^"(e*^^5(e, ■)){x)d^,' 

J—TT 

where is as in (1.5), 

(2.2) 5(e,a:):=E^(^ + 2i^)^''"'"' 

periodic, denotes the Bloch transform of g, g{^) := ^ j'^e~^^^ g{x)dx the Fourier transform, and 

(2.3) g{x)= r e'^-g{C,x)d^ 

J —TV 

the inverse Bloch transform, or Bloch representation oi g & L^. 

^Scc [BJNRZl, BJNRZ2, BJNRZ3, JZN] for related analyses in such general settings. 
^In other words, {e*'^g){^,x) = (e*^«g(5, ■))(»), a consequence of (1.6). 
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Note that, in view of the inverse Bloch transform formula (2.3), the generahzed Hausdorff- Young 
inequahty ||ti||LP(M) < ||^i||L9(^,LP([o,i])) for g < 2 < p and | + | = 1 [JZl], yields for any 1-periodic 
functions /(Cr), C G [-vr,7r], 

(2.4) II r e^«7(C, ■)daLPiR) < ||/||L«(«,Ln[o,i])) for g < 2 < p and - + - = 1.^ 

J —TV P 1 

This convenient formulation is the one by which we will obtain all of our linear estimates. 

Note by (Dl) that there exists a simple zero eigenfunction u' of Lq, which by standard perturba- 
tion results [K] thus bifurcates to an eigenfunction ^(^, •), with associated left eigenfunction ^(^, •) 
and eigenvalue 

(2.5) x{0 = ai^-de + O{\e), 

where a and d are real and d > by assumption (D2) and complex symmetry, A(^) = each 
of (f), 4>, A analytic in ^ and defined for |^| sufEciently small. 

3. Basic linear estimates 
Loosely following [JZl] decompose the solution operator as 

(3.1) S{t) = SP{t) + S{t), SP{t) = u' sP{t), 
with 

(3.2) {sPit)g)ix) = r e^«-a(Oe^(«)*(^(C, ■))LH[o,i])d^, 

J —TV 

and 

iS{t)g){x) := r e^«-(l-a(0)(e^«*5(6)(^)rfe+ T e^«^a(0(e^«*n(05(0)(^)rf^ 

(3 3) J —"rv J —TT 

+ r e*«-«(Oe^(«)^(<^(C,a;) - 0(0, x))(<^(0, 5(0)L2([o,i])dC, 

J —TV 

where a is a smooth cutoff function supported on ^ sufficiently small,^ 
(3.4) nf(O:=0(O('^(O,-)L^([o,i]) 

defined for sufficiently small ^ denotes the eigenprojection onto the eigenspace Range{(^(^)} bifur- 
cating from Range{'u'(x)} at = 0, 4> the associated left eigenfunction, and n := Id — W. 

We begin by reproving the following estimates established (in a slightly different form) in [JZl], 
describing linear behavior under a localized perturbation g G L"'^(IR) n L^(M). 

Proposition 3.1 ([JZl]). Under assumptions (H1)-(H2) and (D1)-(D2), for all t > 0, 2 < p < oo, 
(3 5) d'd^s^im^a < min I + i)"^^'"'''^^"'^ IbllLi(E), 

for Q < n < K + 1, and for some rj > and < / + 2m, n < K + 1, 

(1 + t)"2(^"^/P)~2 ||5f||^i(]j^^j^,+2m+l(IR), 



(3.6) dldrS{t)d:g 



LP(R)-'^'''\e-^'\\d:^g\\Hi+2^+im + {l + t)-^^''^^^^^^ 



Estimates (3.5)-(3.6) were established in [JZl] by first passing to a Green function formulation. 
Here, we give an alternative proof within the Bloch formulation that will be useful for what follows. 



'''Here, and elsewhere, we are adopting the notation ||/||L<!{^,iP([o,i])) := f^^ OlliPdo i])*^?) 
^That is, a = for \^\ > 2^o and a = 1 for \^\ < ^o, sufficiently small. 
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Proof, (i) (Proof of {3.5) (1)). First, expand 



where (^j(^) denotes the jth Fourier coefficient in the Fourier expansion of periodic function ^(^, •), 
and z* = z denotes complex conjugate. 

By Hausdorff-Young's inequahty, |^| < ||5'||li(]r), and by (D2), |e'^'^^)*a^/^(^)| < e~'^^ ^, rj > Q, 
while by by Cauchy-Schwarz' inequality, 



3 V ^ ^ 

Combining these facts, and applying (2.4), we obtain for l/q + l/p = 1 

\\s^it)9\\LP{^;dx) ^ l|e~''^^*||L9{[-7r,7r],(iO ^Up ||^|ki([0,l];dx) I|5'IIl1{M) ^ (1 + *)"^^^~^^^^ II^IIlI (K) ' 

yielding the result for I = m = n = 0. Estimates for general l,m,n > go similarly, passing 
derivatives onto 4> in the inner product using integration by parts and noting that d^. and d^^ 
derivatives bring down harmless bounded factors {i^Y and A(^)"*. 

(ii) (Proof of (3.5) (2)). The second estimate on follows similarly, but substituting for the 
estimate of term /^^ e*^^a(^)e'^*^^)*((^(^), ^(^))^2([o,i])t^^ the slightly simpler estimate 
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e'^''a{Oe^^^^\4>,9)LH[o,i])d^\\Lp(R;dx) < Il«(0e "^^HO, diO) lH[o,i])\\\ Li {i,Lp ([0,1])) 



Lmo,l]-dx)\\Li([-Tr,n];d^) 



< l|e ''^ ^\\L^i{[-7T,7T];d0\\9\\L-^{R;dx) 



where 1/r + 1/s = 1 and qs = 2, so that rq = is oo for p = = 2 and 2 for g = 1, p = oo. This 
verifies the result for Z = m = n = 0; estimates for general l,m,n > go similarly. 

(ii) (Proof of {3.6)(1)). By (D2) and Priiss' Theorem [Pr], we have 

|e^«*(l - a{O)\H'+^[0,l];dx)^H'+'L{[0,l];dx)^ l«(Oe^^*n|^i+i([o,l];(ixHif'+i([0,l];da;) ~ e"''*, ?7 > 0, 

whence, by Sobelev embedding, 

|5^e^«*(l - a{S,))\Hi+i([o^i]-dx)^LP{[0,l]:dx)i \9iaiO^^^*^\H'+^{[0,l];dx)^LP{[0,l];dx) ^ e"''* 

for 2 < p < oo; see [JZl] for details.^ The W^'P(R) norms of the first two terms of (3.3), by (2.4) 
and Parseval's identity, ||5||l2(^^j:/!+iqo,i])) ~ ll5'lli?'+i(M)/'' thus bounded by Ce~^^\\g\\ffi+i(^^y 
The W'^'P{W) norm of the third term may be bounded similarly as in the estimation of above. 



^Here, we are using the fact (see for example the Evans function analysis of [G]) that H''~^^{'M.) and L^(K) spectra 
coincide. 

^'^More precisely, 

llflllHi+lCM) = llffllz,2(5,L2(10,l])) + ll5lli2(£_ij!+l([o,l])) •= llsllz,2(J,Z,2([o,l])) + \\{dx + «0'^^ff(OllL2(5_i2([o,i]))- 
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noting that the factor {4>{$,) —(f)(0)) ~ |^| introduces an additional factor of (1 +t) decay. This 
estabhshes the result for m = n = 0; other cases go similarly, noting that dte^^^fl = L^e^^Il, with 
a second-order operator, so that we may essentially trade one i-derivative for two x-derivatives. 

(ii) (Proof of (3.6) (2)). The second estimate on S follows similarly, but substituting for the 
estimate of term e^«^a(Oe^(^)^(^(C, x) - ^!>(0, a;))(^(0, 5(0)l2([o,i])^^C the estimate 

II r e^«-a(Oe^(«)*(0(C,x) - <^(0, a;))(<^(0, 5(0)^^ellL.(.) < ll«(Oe-^'^'*l€ll(<^(0,ff(0)lllL.(^,L.(.)) 

J —TV 

<llie|e-''«'*||KOIlL^(.)llL.(0 
< llie|e-''«'*||L^.(oll5llL^(.) 

<{l+t)-hW^-yP)-h\\g\\,,, 

where 1/r + 1/s = 1 and qs = 2, so that rq = is oo for p = g = 2 and 2 for g = 1, p = oo. □ 
Remark 3.2. In [JZl], by contrast, sP{t) is estimated through its integral kernel 

e{x,t;y) := = T a(Oe^(«)*e*«^0(e, y)^^, 

J— IT 

yielding the slightly sharper estimate < snpy \\e{- , t; y)\\ < ||0|koc(,,^)(l+t)-V2(i-i/p). 

For our purposes, this makes no difference, and, as we shall see in Section 4, there can be an advan- 
tage in maintaining the separation into distinct frequencies afforded by the Bloch representation. 



4. Linear behavior for modulational data 

Next, we consider behavior of (2.1) when applied to modulational data g = h^v! , the linearized 
version of a nonlinear modulational perturbation u{x + ho{x)) — u{x) ^ ho{x)u'{x). The following 
estimates, obtained by frequency-domain rather than spatial-domain (Green function) techniques 
as in [JZl, JZ2, JZN], together with the associated modified nonlinear iteration scheme of Section 
7, below, represent the main new technical contributions of this paper. 

Proposition 4.1. Under assumptions (H1)-(H2) and (D1)-(D2), for allt>0, 2 < p < oo, 

(4.1) \\didrsnt){hou')hr>^^^ < (1 +t)-|(l-l/f)+|-^||5./l0||Li(R), 

for I + m > 1 or else I = m = and p = oo, and, for < I + 2m < K + 1, 

(4.2) \\dldrS{t){hou')h.^^^ < {l + t)-"^^'-'/P^\\d,ho\\mm)nHi+^m+HR), 
and when t < 1 

\\di:di^{SP{t) -ld){hou')\\Lp(^) < ||9^/iollLi(R)n//'+2'"+i(]R)) 
\\dicd^{sP{t){hou') - ho)\\Lp(^^) < ||9^/io||Li(M)nL2(M)- 



7 



Proof, (i) (Proof of (4.1)j, I + m > 1. We treat the case Z = 1, m = 0; other cases go similarly. 
First, re-express 

d,isPit){u'ho))ix) = r iCa{Oe^^^^'e'i^^u',ho)m[o,i]){Od^ 

J —n 

= Y. r i^o^iOe'^^^'e'^^4>{^, y)u\y), e^'^^^) L^^^o,i]Mi + 2j7r)dC 

(4.4) 

where 4'{S,)u'j denotes the jth Fourier coefficient in the Fourier expansion of periodic function 

m^'- By (2.4), we thus get ||9x(sf (t)(^i'/io))||Lp(E) < C{1 + t)-"2i^-^M\\d^ho\\Li(R)j:^ 
yielding the result together with the Cauchy-Schwarz estimate 

I d>(£)u' ■ I — — — 

3 V i J 

(ii) (Proof of (4.1)j, Z = m = 0. This follows by an explicit error function decomposition 
as in [JZ2], p. 18, putting a ^ factor on h^u' and a 1/^ factor on e^*^^)* and proceeding as in 
the cited estimate. More precisely, for the low-frequency part of Ziq, express the principal part 

+0O , (aiS,-de_'^)t 



/+0O (ml;-dt,^)t 
^i^x^ ^ errfn((x - atf/Vt) 
-oo K 

mmus 

^./|C|>7r e^^"^ 7^1/2 estimating the L°°(M) norm of convolutions with the former by the 

Triangle Inequality, and with the latter by Hausdorff Young's inequality; for details, see [JZ2].^^ 
(in) (Proof of (4.2) ). Likewise, this follows by re-expressing 

S{t){hou'){x) = r e'^^l-a{0){e'^^%{Ou'){x)dC+ f e'^-a{i){e'^^'m)h^{i)u'){x)di 

J —n J —IT 

+ r e^«^a(Oe^«)* (</)(?, x) - mx)){mMi)u') L\[o,i])di 



(4.5) f^J^ sft + a^j) 



^^This estimate is not needed for the nonlinear iteration, so we do not emphasize it here. We note that the 
resulting bounds (1.8) may be recovered alternatively by the modulation decomposition of [JNRZl] without case 
I = 0. 



then estimating as before, where we are using (1 — a(^)) < ^, n(0^' = 0{^) , — = 0(C), 

respectively, to bound the key terms ll[+2nj) ' ^'^fii+L]) ^ ' 1(|+2tjO ^ ^ ^' ^"^^ ^^^"^^ 
the Cauchy-Schwarz inequahty 

dxhoi^ + 2j7r 



to bound the sum over j ^ 0. This estabhshes the case I = m = 0; other cases go similarly. 

(iv) (Proof of (4.3) ). Expanding SP{t) - Id = {SP{t) - 5^(0)) - 5(0) = tdtSP{s{t)) - 5(0) for 
some < s{t) < t, we obtain the first inequality by combining (4.1) and (4.2). Likewise, the second 
ineqality follows by expanding sP{t){hou') — ho = {sP{t) — sP{0)){hou') + {sP{0){hou') — ho) and 
applying (4.1) together with 

\\sP{0){hou') - ho\\Lp{R) 

< CsupgEj |g+2j^| |Px/to|Ui(R) + <^SUp^Ej^O |g+2j7r| ^ <^|Px/lO || Li (ffi)nL2(ffi) 

This establishes the case Z = m = 0; other cases go similarly. □ 

Remark 4.2. 1. If we split /to into high-frequency and low-frequency parts, then the contribution 
of the high-frequency part decays faster by factor (1 + t)~^l'^ in all estimates (4.1)-(4.3). 

2. In the estimate of s^, it is easy to see that the same bounds hold if u' is replaced by any 
periodic / G iJpgj.([0, 1]). However, in the estimate of S, replacing v! by periodic / G i?pgr([0, 1]) 
introduces an exponentially decaying error in / — reflecting the dynamics of the problem on a 
periodic domain. To extend our results to the sum of a localized perturbaton and a perturbation 
asymptotic as x — >■ ±oo to any two fixed periodic waves, not necessarily modulations, requires only 
the estimation of this latter error term as a time-exponentially function from L°°(R) — > L°°(]R), a 
semigroup/Fourier multiplier problem of concrete technical nature. This would be an interesting 
direction for further investigation. 

3. The estimate (4.3) is the key to handling the "initial layer problem" in our later nonlinear 
iteration, allowing us to essentially prescribe initial data for ^ as is convenient for the analysis; see 
Remark 7.2 below for more information. 

5. Nonlinear perturbation equations 

Essentially following [JZl], for u{x,t) satisfying ki,ut = k^Uxx + f{u) + ki,cux and tpi^-it) to be 
determined, set 

(5.1) u{x,t) = u{x — 'ip{x,t),t) and v{x,t) = u{x,t) — u{x). 
Lemma 5.1 ([JZl]). The nonlinear residual v defined in (5.1) satisfies 

(5.2) k^{dt-L){v + ijUx) = k^M, kj\f =Q + nx + {kA + kldl)S + T, 
where 

(5.3) Q:= f{v + u)- f{u)-df{u)v 



(5.4) Tl := -Kvipt - k^vtl^xx + Ki^x + 



Xj -. I ■ 



(5.5) 



S := vipx 
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and 

(5.6) T:=-{fiv + u)-f{u))^^. 

Proof. (See also [JZl].^^) Prom definition (5.1) and the fact that u satisfies (1.2), we obtain 

(5.7) A;*(l - tl)x)ut + K{-c + tlJt)ux = kl {^ ^^^ + (1 " '^x)f{u) 
hence, subtracting the profile equation —k^cux = k^Uxx + f{u) for u, 

(5.8) M. - KLv + - k,i,.n + Km. = a - i>.f(n) + kl {j^u. 

or 

(5.9) fc^(9t - L)v + fc^(Vn:,)t = Q + + {kA + ^^9^)5 - V^/(n) + kl {tpxUx)x , 
and we finish the proof using Lux = and the profile equation for u to obtain 

(5.10) k^L{tpUx) = k-^c^xUx + kl'^xUxx + kl {'4^xUx)x = -tpxf{u) + kl {'^xUx)x ■ 



□ 



6. Nonlinear damping estimate 

Proposition 6.1 ([JZl]). Assuming {HI) - {H2), let v(-,0) G {M) (for v as in (5.1)) and 
suppose that for some T > 0, the H^{M.) norm of v{t) and ipt{t) o^c^ the H^^^{M.) norm of ipx{t) 
remain bounded by a sufficiently small constant for allQ <t <T. Then there are positive constant 
9 and C, independent of T, such that, for all <t <T, 

(6.1) ||i;(t)||^K(K) < Ce-''\\vml^^^^ + C re-^(*-^) {\\vis)\\l,^^^ + ||(Vt, V'.)(^)||^k(k)) ds. 

J 

Proof. Take for writing simplicity A;* = 1. Rewriting (5.9) as 

(1 - '^^)vt - Vxx - cvx = df{u)v + Q- {Ux + Vx)ll^t 

(6.2) _ / _ ^2 

+ {{Ux + Vx)iIJx)x + [{Ux + Vx)- ^ - f{u + V)lpx, 

\ 1 — W-r^ X 



taking the inner product against X^jLo ^"T^^^^' integrating by parts, and rearranging, we obtain 



(t)<-^l5f+Vi)||i.(K) + c(|| 

^(*)ll?f^(R) + ll(V't,^x)(t)|llfif(K)) , 



for some > 0, so long as ||(f , V't) V'x) V'xx)(i)||H^f+i(R) remains sufficiently small. Sobolev interpo- 
lation Ibll^Km) < C""-^ 119^+^511^2(5;) + (51151112^) gives, then, for C > sufficiently large. 



d 

lull 



^^n-.l|^(R)(*) < -^ll^'(*)ll|i^(E) +C (ll^Wlli^W + ll(V't,V'.)(i)llHK(K)) , 

from which (6.1) follows by Gronwall's inequality. See [JZl] for further details. □ 



Note that we have here followed a different convention than in [JZl], reversing the sign of ip in (5.1) in agreement 
with formal asymptotics of [Se, DSSS]. The change — >■ — recovers the formulae of [JZl]. 
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7. Nonlinear iteration scheme 



The key idea is, similarly as in the localized case treated in [JZl], starting with 

{dt - L){v + il)u') = J\f, v\t=o = do, tl'\t=o = ho, 

where do := wo(- - ho{-)) - ue L^{R) n H^{R), d^ho G L^{R) n H^(R), to choose ^ to cancel 
contributions, as 



(7.1 



^(i) = sP{t){hou + do) + f sP{t - s)Af{s)ds 

Jo 



- (1 - x{t)) (^sP{t){do + hov!) -ho + sP{t - s)M{s)ds 
where x(t) is a smooth cutoff that is zero for t < 1/2 and one for i > 1, leaving the system 

v{t) = S{t){do + hou')+ j S{t-s)M{s)ds 

(7-2) . 

+ (1 - x{t)) \ySP{t)do + {SP{t) - ld)hou' + j SP(t - s)Af{s)d 

We may extract from (7.1)-(7.2) a closed system in {v,iljx,'(pt) (and some of their derivatives), and 
then recover ijj through the slaved equation (7.1). 

Remark 7.1. At first sight, we have accomplished nothing by introducing a ^'-dependent change of 
variable and choosing ■0(0) = since we have replaced the nonlocalized perturbation uo{x) — u{x) 
used in the previous ho = setup of [JZl], by a different nonlocalized perturbation do + hou'. 
However, what we really did was replace the asymptotic states u{x + c±) — u{x) by their linear 
approximates c±u'{x), which removes the key difficulty of higher order remainders in the Taylor 
expansion of nonlinear modulations. 

Remark 7.2. Notice that modulational data u'ho enters in (7.2) only through operators S{t) and (1— 
x{t)){SP{t) — Id) for which we have Gaussian decay in LP(R) with respect to ||c^x^o||Li(R)n//i{R)) 
hence the error incurred by defining iJj by (7.1) instead of the value iplt) = s^{t){hou' + do) + 
Jq sP{t — s)J\f{s)ds exactly canceling terms is harmless to our analysis. The choice of (7.1) 
reflects our need to accomodate the incompatibility between the initial value tp\t=o = ho prescribed 
by the spatially-asymptotic behavior of the initial perturbation and the function encoding time- 
asymptotic behavior of the perturbed solution; that is, it is a device to avoid having to resolve an 
initial layer near t = 0.^^ Whether this initial layer is an artifact of our analysis or reflects some 
aspect of short-time behavior is unclear; as the estimates show, this is below our level of resolution. 



8. Nonlinear stability 

With these preparations, the proof of stability now goes essentially as in the localized conservative 
case treated in [JZ2], using the new linear modulation bounds to estimate the new linear term 
coming from data Hqu' in (7.1) and (7.2). As noted in [JZl], from differential equation (5.2) 
together with integral equation (7.2) (7.1), we readily obtain short-time existence and continuity 
with respect to t of solution {v,'tpt,'<px) € iJ^(M) by a standard contraction-mapping argument 



For simplicity we set Z = m = in this discussion. 
■'■^In the case ho = 0, this essentially reduces to the simpler device used in [JZl] to treat the localized case, of 
substituting x{t)s''{t) for s^{t). However, the latter is clearly too crude to treat the present case. 
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treating the linear df{u)v term of the lefthand side along with Q,Tl,S, T, ipu' terms of the righthand 
side as sources in the heat equation. Associated with this solution define so long as it is finite, 

(8.1) C{t)-= sup \\{v,tpt,'(px)is)\\HKm{l + s)^/'^. 

^ ' 0<s<t 

Lemma 8.1. For all t > for which (^{t) is finite and sufficiently small, some C > 0, and 
Eq := 11(^0, 5a;^o)||Li(M)nif^(iR) Sufficiently small, 

(8.2) at)<c{E^+at?)- 

Proof. By (5.3)-(5.6) and corresponding bounds on the derivatives together with definition (8.1), 
and using (6.2) to bound vt, 

(8.3) ||Ar(t)|Ui(M)ni/i(R) < ll(^,V't,V'x)(i)ll|3(K) < CC(t)'(l + t)-^ 

so long as C,{t) remains small. Applying the bounds (3.5) (l)-(3.6) (1) and (4.1)-(4.3) of Propositions 
3.1 and 4.1 to system (7.1)- (7.2), we obtain for any 2 < p < oo 

, , ||^;(i)||Lp(R) <C(l + t)-5a-Vp)£;„ + C'C(t)' / (l + «-s)-5(i-Vf)-i(i + s)-|ds 

(8.4) Jo 

<Cp(£;o + C(i)')(i + *)"^^'"'/^^ 

and 

||(V't,V.)(i)||w^™(M)<C(l + t)-^£:o + CC(t)2 /■\l + t_5)-^(i-i/f)-i/2(i + ,)-irf, 

( o.o 1 ^0 

<C^{E^^C{tf){\^t)-\^^-^l^\ 

Estimate (8.5) yields in particular that ||(V't) V'a;)(*)||//^+i(R) is small, verifying the hypotheses of 
Proposition 6.1. From (6.1) and (8.4)-(8.5), we thus obtain 



<C(i?o + C(t)')(l + i) 



Combining this with (8.5), p = 2, rearranging, and recalling definition (8.1), we obtain the result. 

□ 

Proof of Theorem 1.1. By short-time existence theory, ||(f, V'i; V'x)(i)llj/^(R) is continuous so 
long as it remains small, hence C remains continuous so long as it remains small. By (8.1), therefore, 
it follows by continuous induction that C,{t) < 2CEq for t > 0, if Eq < 1/4C, yielding by (8.1) the 
result (1.7) for p = 2. Applying (8.4)-(8.5), we obtain (1.7) for 2 < p < for any < oo, with 
uniform constant C. Taking p* > 4 and estimating 

(8.6) \\{Q,TZ,S,Tm\\LHR) < \\{v,iPt,4^.,4^..m\\%^^) < CEo{l + t)-i 

in place of the weaker (8.3), then applying (3.5) (ii) in place of (3.5) (i), wc obtain^^ 

\m, ^x)(t)iivK-+i-w ^ ^-(1 + tr'^Eo + catf j\i + 1 - ,)-^(i/2-i/P)-i/2(i + ,)-!^, 

<C(i?o + C(i)')(l+i)-^^'-'/^\ 



^"'Compare to the argument of Lemma 4.2, [JZ2], regarding localized perturbations in the conservative case. 
"'^''We bound the dtS contribution according to 

/* s^t - s)dtS{s)ds = - f dt[s'']{t - s)S{s)ds + s''{0)S{t) - s^(t)>S(0). 
Jo Jo 
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for 2 < p < oo. Likewise, using (8.6) together with bound 

ii(Q,r)(i)iki(M) + wmwH^R) + ii'5(i)iii/3(M) < c(i)'(i +t)-^ 

obtained in the course of proving (8.3), we may use (3.6) (ii) rather than (3.6) (i) to get 

MmLPiR) < C{1 + t)-^2('-'MEo + Catf fe-'^^'-'Xl + s)-"^ds 

Jo 

+catf r(i+t-s)-5(V2-i/f)-^(i+,)-!ds 

Jo 



<c{Eo + citf)ii + ty 



and achieve the proof of (1.7) for 2 < p < oo. 

Estimate (1.8) then follows through (7.1) using (3.5) (i), by 

(8.9) Wmh^m < CEo + cat? / (1 + i - s)-^(l + s)-'2ds < C{Eo + Citf), 

Jo 

yielding nonlinear stability by the fact that 

(8.10) u{x — ip{x, t), t) — u{x — ip{x, t)) = v{x, t) + u{x) — u{x — ip{x, t)), 

so that \u{t) — u\ is controlled by the sum of \v{t)\ and \u — u{- — 'tp{-,t))\ < \ip{t)\ sup 

Acknowledgement. Thanks to Bjorn Sandstede for pointing out the results of [SSSU]. 



□ 
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